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FILE  TRAHSi^  FUNCTIONS 


The  purpose  of  this  memo  is  to  present  a  detailed  derivation  of  the 
transfer  function  of  a  critical  pile  and  to  relate  this  steady-state  description 
and  the  transient  response  to  a  step  fxinction  change  of  reactivity.  Reference 
to  G.E,  report  JIO-1  -  Revised  is  implied  throughout  this  discussion. 


The  transfer  function  of  a  critical  pile  can 
differential  eqiiations  of  the  pile.  The  reactivity  of 
he  controlled  hy  the  delayed  neutrons  so  that  keff 
For  these  assumptions  elementary  pile  theory  gives  - 
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<1  °l/S  ^1  \ 

— - n  -  Aj  C^/B 


he  determined  from  the 
the  pile  is  assiuned  to 
always  less  them  1  +  B 

(1) 

(2) 


where ;  f^ 
n 

/ 

B 

Ci 

Ai 

t 

D 

S^eff 

keff 


=  Bi/B 

s  neutron  flux 
s  mean  life  of  a  neutron 
s  ^  of  group-i  fission  neutrons 
=  ^  of  neutrons  emitted  which  are  delayed 
r  concentration  of  fission  fragments  which  emit 
group-i  delayed  neutrons 
s  decay  constant  of  group-i  neutrons 
s  time  (seconds) 

=  reactivity  in  dollars  - 

=  keff  -  1 

a  ratio  of  fission  neutrons  produced  to  those  lost  for  an  infinite  pile 


^^eff 
B 


The  discussion  thus  far  is  somewhat  misleading.  We  have  been  con¬ 


sidering  a  pile  on  the  basis  of  elementary  pile  theory,  when  actually  we  are 


going  to  use  an  electronic  pile  simulator  whose  correspondence  to  a  pile  de¬ 
pends  on  the  accuracy  of  the  assumptions  of  elementary  pile  theory.  To  avoid 
considering  these  assumptions  we  will  he  more  exact  and  say  that  we  desire  to 
determine  the  transfer  function  of  the  electronic  pile  simulator.  The  simulator 
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obeys  the  same  differential  equations.  To  determine  the  transfer  function  consider 
the  simulator  as  a  four-terminal  network.  The  input  function  is  the  reactivity^ 

D  =  5  (^)  the  output  is  the  neutron  flux.  For  a  sinusoidal  function  of  time 

input  the  transfer  function  E(j(e)  is  the  ratio  of  the  Laplace  transforms  of  the 
output  and  input.  It  consists  ot  K  -  a  constant,  or  frequency  invariant  portion, 
and  &(s)  s  the  frequency  variant  portion,  or  the  portion  dependent  on  time 
derivatives.  In  general,  we  define  fCs)  as  the  Laplace  transform  of  f(t)  by 
this  integral:  ^ 

f(s)  s  (e"®*  f(t)  dt 


S(t) 

n(t)  Sfe) 

F&(s) 

n'(s) 

^(s)  K6(s)  =  n'Cs) 


Let  the  reactivity  nave  small  sinusoidal  oscillations  of  amplitude  ( 

D  =  S(t)  =  ee  (3) 

Since  the  pile  is  critical,  the  neutron  flux  will  vary  sinusoidally  n*  about 
a  steady  level,  Ug.  The  same  is  true  for  the  concentrations  of  fission  fragments. 


Thus: 

n(t)  =  no  n'(t)  (if) 

Cj^Ct)  =  Cj^o  /  ®i'  (5) 


By  definition  the  transfer  function  relates  sinusoidal  variations,  so  that 

r'-UI 

K&(s)  = 

To  obtain  this  result  we  eliminate  from  equations  (l)  and  (2), 
take  the  Laplace  transform  and  substitute  to  obtain  an  expression  for  n'. 

From  (2)  ^i°i  _  £i 

B  T  dt 


(6) 


(7) 
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Substituting  in  (l)  and  using  (3) 


I  (g_  i)n  /Xfi  n-Z/dCC^ 


By  (4)5(5): 


d. 

B  dt 


=  nQS  “i 

Y 

but  n  =  ■^  f  j^n 


i  dt 


It  is  convenient  to  neglect  n'(t)  S(t).  Bor  power  levels  greater  than  one 
quarter  power  Ug ^ (t))$>  n'  (t)  s  (t),  so  that  this  assumption  vfill  not  introduce 
appreciable  error.  Taking  the  Laplace  transform  of  (9)  gives 


s  n' (s)  =  no  S (s) 


C.  (£ 

i  B  1 


Rewriting  eq.  (2)  in  terms  of  n*(t)  and  Cj^*(t)  we  have: 

d  Cj^/B  f^  /  fi  \  • 

=  7  “0  i  4“  -  -3- 

But  the  steady-state  condition  dCj^/B  _  q  shows  that 

dt 

^i  “0  X  ^Iq 

taking  the  Laplace  transform  and  then  solving  for 

B  ^  B 

C^  » (s)  _  f^  n*(3) 

B  ^  (s  /  ^i) 

Substituting  (13)  in  (lO) 


I^LL3l  -  ^0 _ —  =  ajL^r  — 

S(s)  SiiT.  fj  SUT  A 


- - 

1  iT  hjJ. 

^  ^r(s  /M) 


(15) 
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SJhis  can  be  rewritten: 


=  n„B 


(8  /Xi)(s  /Ap)...(s 


si^.  =  ^  _  (17) 

(s)  Sq  7^  ais  /  ^2®^  ^••*  /  f  ••• 

=  52®  ^,.^.®.  -  (18) 

s(s  /  s-j^)(s  /  S2)»*<'(s  /  Sj^)  ••• 

This  is  the  general  form  of  the  transfer  function.  All  the  Aj^'s  and  s^'s  are 
positive.  That  Sj^  is  positive  follows  from  the  fact  that  all  the  a^^'s  being 

positive  ezcludes  positive  roots.  This  just  says  that  the  pile  transfer  function 
can  be  handled  by  the  standard  techniques  of  servo  theory.  (This  was  implicit 
in  the  original  assumptions). 

This  general  result  may  now  be  applied  to  a  simulator  having  five 


groups  of  delayed  neutrons, 
B<J.f-life 


0.05 

0.43 

1.52 

4.51 

22.0 


Assiuae 


Decay  constant 

X i  ■  sec“^ 

*1/2 

XI  =  14 
X2  r  1.6l 
X3  =  0.456 
X4  =  0.154 
X5  =  0.0315 

^  -  5  at  10“^  sec 


Fraction  of  Mission  Neutrons 


=  0.00029 

B2  =  0.00084 
B3  s  0.0024 
B4  =  0.0021 

B3  =  0.0017 


Eq.  (17)  is  now 

n*(s)  _  UpB  (s  /Xi)(s  f  X2)(5  /X3)(s  f  X4)(s  /X5) 

(s)  jt  %  s5  /  aj[j,s^  j-  /  ^2®^  ^  ^  ®'0 

The  a. 's  can  be  written  in  terms  of  Bj^  and  Yby  comparing  them  with 


the  corresponding  coefficients  of  the  denominator  expressed  in  eq.  (16). 
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The  results  are  shown  in  Appendix  I.  Now  the  roots  of  this  fifth  degree  eq.uation 
must  be  determined. 

D(s)  B  s5  /  163  s^  /  2270  s3  /  3750  s^  j-  1200  s  /  73.2  =  0  ^20) 

=  (s  /  s^)(s  f  S2)(s  /  ^-^{5  ^  e^)(B  -f-  s^)  (21) 

The  results  are:  S]^  =  148  , 

sg  r  13.3 

53  -  1.45 

54  =  0.32 

S3  r  0.08 

These  results  can  be  used  to  plot  the  Bode  diagram  of  figure  1  by  replacing  s 
by  jw. 

The  Bode  diagram  (Pigure  l)  shows  the  asymptotes  of  the  steady-state 
amplitude  gain  of  the  pile  simulator  for  a  sinusoidal  input.  Eeplacing  s  by  jw 
is  an  obvious  operation,  since  s  may  be  considered  equivalent  to  the  time 
derivative.  Taking  the  time  derivative  of  a  sinusoid,  e:^  jwt,  is  just  multiplicar- 
tlon  by  jw.  The  values  of  Sj^  and  are  the  break  points  of  the  asymptotes  of  the 


function  whose  analytical  expression  is; 
Lm|K&(jw)|  B  Lmn^  /  Lm  ^  -  Ian 

/  Im  j  jw  /  l4|  -*  Ian 

/  Ian  I  jw  /  1.6l[  -  Ian 

/  Ian  I  jw  /  0.456|  -  Ian 

^  Ian  I  jw  0.l54j  •  Ian 

/  Lm  I  jw  0.0315|  -  Ian 


Transient  Response 


jw  ^  l48j 

>  /  13.3 

jw  /  1.45 1 

jw  /  0.32 

jw  /  0.08| 


The  response  of  a  steady-state  pile  to  a  step  function  change  of 
reactivity,  S  ,  can  be  expressed  as  a  sum  of  exponential  terms.  Usually 


N  B  N. 


^igPi 
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we  are  interested  in  six  groups  of  delayed  neutrons  giving  seven  terms  in  the 
eacpaasion.  For  each  amplitude,  S  ,  there  are  seven  values  of  p.  Each  p  has  an 
A  associated  with  it.  Eigure  2  shows  the  calculated  values  of  Pj^  and  A^  for 

values  of  reactivity,  S  ,  from  ~  $2.00  to  -f  $1.00  (-1.5^  to  -f-  0.75^).  These 
curves  also  are  plotted  for  various  values  of  the  mean  life  of  a  neutron,  J. 

These  step-function  response  graphs  can  be  related  to  the  Bode 
diagram  of  steady-state  sinusoidal  response. 

The  transfer  function  of  the  simulator  KQ(s)  depends  only  on  the 
simulator's  physical  characteristics.  On  the  other  hand  the  constants  obtained 
from  the  step-function  response  graphs  depend  on  the  amplitude  of  the  input. 

If  we  apply  a  unit  impulse  function  (unit  step-function  minus  another  unit  step- 
function)  we  can  use  the  step-response  graphs  to  determine  the  output.  The  values 
of  Aj^  and  pj^  are  those  corresponding  to  iii  a  sense  do  not  depend  on 

reactivity  S  ,  but  only  on  the  characteristics  of  the  simulator.  Now  by  previous 
definitions  - 

n(s)  s  S(s)K&(s)  (2h) 

If  S(t)  is  an  impulse  function  5(®)  =  1 


n(3) 


(t)=0 


The  step-function  response 

6 

n(t)  a  Uq  Aj^  Pi^ 
i=l 


has  the  Laplace  transform 

/  N  ^i 

n(e)anQ^  — 

i-1 


Pi 


(25) 


(26) 


(27) 


subject  to  certain  mathematical  restrictions  (the  physical  equivalent  of  not 
allowing  the  pile  to  blow  up).  Eor  the  case  of  6  g^o^^ps  of  delayed  neutrons 
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the  transfer  function 

E}(s)  _  ISa  (s  ^\){s  /\)(s  /Kq)(s 


J  s(s  /  S2^)(s  /  S2)(s  /  S2)(s  f  sj^)(s  /  s^)(s  /  Sg) 


Equation  (25)  states  that  this  equals 

1  „  =  *^0  S  — 


(t)=0  “  °  i_i  s  -  pj, 

^  /  hps^  /  hi  s  hf 

°  (s-p^)(s-p2)(s-p^)(s-pj^)(s-p^)(s-pg)(s-p^) 


There  is  a  correspondence  between  the  s^  's  and  the  Pi  ’s. 


Pi  =  o 

P2  =  -0.O15 


s6  =  0.015 


p^  =  -0.08 

S3  =  0.08 

P4  =  -0.33 

»4  =  0*32 

-  -1.45 

S3  =  1.45 

P6  =  -13.3 

S2  =  13.3 

P7=  — 

Sj^  s  148 

The  value  of  sg  was  obtained  from  the  Bode  diagram  of  G.E.  report  JIO-1, 
The  values  of  P]^  to  pg  and  s^  to  S2  do  not  change  appreciably  for  different  values 
of  neutron  life  so  that  the  above  comparison  is  valid  even  thou^  the  p^'s 
were  computed  for^s  10“^  and  the  sj^'s  iox -  5  2:  10“^.  The  value  of  s^  (and  pr,j 
varies  considerably  for  different  values  of  but  for  a  given  ^ ^  Sj^  is  the 


negative  of  p^.  This  correspondence  between  steady-state  and  transient  response 
gives  a  convenient  check  and  can  be  used  to  correlate  results. 
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<^4=  ^ 


+-|^(\'‘'\ +^3  +  ^y)  +^(x,  +  X^  +  \+\) 


ZHS  3  is- 


+^f  \3  ■"  •"  K  ^^=  ■"  M 

+  |l(X„  +  \  *K*\i*\s^^  +f 


I 

t' 


\x3>f  ■*■  ^3^5  |^(  A„,A3^  + 

+-^(Vs +\z5 +  -^-A^vs) 

+  ^(^/i3'^'^a5'^\35+A^3^  -f 


^0““  ‘^lZ3‘is'^'^\3V5-'^^^'3H5  ^Ii4y 


where  \jk  ,..  -  \  \  '  " 


